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Abstract 

The current status of optical potentials employed in the prediction of 
thermonuclear reaction rates for astrophysics in the Hauser-Feshbach 
formalism is discussed. Special emphasis is put on a+nucleus potentials. 
Further experimental efforts are motivated. 

1 Introduction 

The investigation of explosive nuclear burning in astrophysical environments 
is a challenge for both theoretical and experimental nuclear physicists. Highly 
unstable nuclei are produced in such processes which again can be targets for 
subsequent reactions. Cross sections and astrophysical reaction rates for a 
large number of nuclei are required to perform complete network calculations 
which take into account all possible reaction links and do not postulate a priori 
simplifications. 

The majority of reactions can be described in the framework of the statisti- 
cal model (compound nucleus mechanism, Hauser-Feshbach approach, HF) [|TJ , 
provided that the level density of the compound nucleus is sufficiently large 
in the contributing energy window 0. In astrophysical applications usually 
different aspects are emphasized than in pure nuclear physics investigations. 
Many of the latter in this long and well established field were focused on 
specific reactions, where all or most "ingredients", like optical potentials for 
particle transmission coefficients, level densities, resonance energies and widths 
of giant resonances to be implemented in predicting El and Ml 7-transitions, 
were deduced from experiments. As long as the statistical model prerequisites 
are met, this will produce highly accurate cross sections. For the majority 
of nuclei in astrophysical applications such information is not available. The 
real challenge is thus not the well-established statistical model, but rather to 
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provide all these necessary ingredients in as reliable a way as possible, also for 
nuclei where none of such information is available. 



2 Transmission Coefficients 

The final quantities entering the expression for the cross section in the statis- 
tical model [|IJ are the averaged transmission coefficients. They do not reflect 
a resonance behavior but rather describe absorption via an imaginary part in 
the (optical) nucleon-nucleus potential ||. In astrophysics, usually reactions 
induced by light projectiles (neutrons, protons, a particles) are most impor- 
tant. Global optical potentials are quite well defined for neutrons and protons. 
It was shown M, |5j that the best fit of s-wave neutron strength functions is 
obtained with the optical potential by ||, based on microscopic infinite nuclear 
matter calculations for a given density, applied with a local density approxima- 
tion. It includes corrections of the imaginary part [7|, §|. A similar description 
is used for protons. Deformed nuclei are treated by an effective spherical po- 
tential of equal volume ||. For a detailed description of the formalism used 
to calculate El and Ml 7-transmission coefficients and how to include width 
fluctuation corrections, see 0, |5| and references therein. 



2.1 a+Nucleus Potentials 



Currently, there are only few global parametrizations for optical a+nucleus 
potentials at astrophysical energies. Most global potentials are of the Saxon- 
Woods form, parametrized at energies above about 70 MeV, e.g. P, |10|. The 



high Coloumb barrier makes a direct experimental approach very difficult at 
low energies. More recently, there were attempts to extend those parametriza- 



tions to energies below 70 MeV |TT 



Early astrophysical statistical model calculations [12, 13 1 made use of sim- 



plified equivalent square well potentials and the black nucleus approximation. 
Improved calculations jnj] employed a phenomenological Saxon-Woods po- 
tential flql , based on extensive data [pf! - This potential is an energy- and 



mass-independent mean potential. However, especially at low energies the 
imaginary part of the potential should be highly energy-dependent. 

Most recent experimental investigations [T7|, |ll| found a systematic mass- 
and energy-dependence and were very successful in describing experimental 
scattering data, as well as bound and quasi-bound states and B(E2) values, 
with folding potentials [|19|]. Motivated by that description (a global para- 
metrization is not given in [18[]), the following global a+nucleus potential is 



proposed. The real part V = \V{ of the optical potential is calculated by a 
double-folding procedure: 



Vf(r) = / / p P (r P ) p T (r T ) v eS (E, p N = p P + p T , s=\f+fp-fr\) d 3 r P d 3 r T , 

(1) 

where pp, px are the nuclear densities of projectile and target, respectively, and 
v e s is the effective nucleon-nucleon interaction taken in the well-established 



DDM3Y parametrization pOfl . The remaining strength parameter A can be 



determined from the systematics |17j of the volume integral per interacting 
particle pair 

Jfl=-rV / V(r)d 3 r . (2) 
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Introducing a mass-dependence in addition to the energy-dependence, one 
can reproduce the behavior derived from the imaginary part by applying the 
dispersion relation ||18||, by the relation 



J, 



R 



' f(A) + 0.67£ c . m . (E c . m . < 26) 

[0.011/(A) - 2.847] E cm - 1.277 / (A) + 626.591 (26 < E cm < 120) , 

(3) 

with f(A) = 311.0132 exp (A T 2/3 ) and E c . m . given in MeV. Using Eqs. | and 
^, the value of A can be set. 

The mass- and energy-dependence of the imaginary part is more crucial 
for the calculation of transmission coefficients. The volume integral of the 



imaginary part can be parametrized according to |2T 



f for E c m < E 

where E is the threshold energy for inelastic channels, Jo is a saturation 
parameter and A the rise parameter. 



As can be seen |18[], the resulting volume integrals are quite diverse for 



different nuclei and a straight-forward mass-dependent parametrization is not 



possible. In Ref. \\L8\ , Jo and A were fitted to experimental data. For a global 
parametrization, one has to include nuclear structure and deformation infor- 
mation into the fit. Having chosen the (probably energy-dependent) geometry 
parameters |l]J (and appropriately modified them for deformed nuclei), the 
depth W of an imaginary Saxon-Woods potential can be determined with 
Eq. £|. At energies E C Tn > 100 MeV, the value of W should reach the satura- 
tion value [I0"|. This determines the parameter J . The threshold E Q = E (p) 



is given by the energy of the first excited state. The increase in the num- 
ber of inelastic (competing) channels is described by A. As a first estimate, 
this can be related to the increase in the total number n = n(p) of levels: 
A = / (dn/dE) = aE^ + (3A 2 ^ + 7 In (dn/dE). The functional dependence is 



found by comparison with experimental values for A ||18|| . Thus, microscopic 
and deformation effects are implicitly included via the level density p = p(E). 
With nuclear level densities p{E) taken from ||, a good fit is obtained to the 



values given in fT5fl , as well as for various (0,7) and (n,a) data P^fl . However, 



it has to be emphasized that the aim is to obtain a reasonable global potential, 
i.e. with a low average deviation from experiment. For specific reactions, such 



a potential might be improved by fine-tuning to a particular nucleus (e.g. |23|, 
but will lose predictive power by that. 



3 Conclusion 

With the recent improvements of the level density treatment P] and the 
a+nucleus potential, the two major remaining uncertainties in the existing 
global astrophysical reaction rate calculations have been attacked. The new 
descriptions will provide more reliable predictions of the low-energy cross sec- 
tions of unstable nuclei. 

Nevertheless, to check and further improve current parametrizations, not 
only of a+nucleus potentials but of all involved quantities, experimental data 
is needed. Especially investigations over a large mass range would prove useful 
to fill in gaps in the knowledge of the nuclear structure of many isotopes 
and to construct more powerful parameter systematics. Such investigations 
should include neutron- and proton-strength functions, as well as radiative 
widths, and charged particle scattering and reaction cross sections for stable 
and unstable isotopes. This information can be used to make future large-scale 
statistical model calculations even more accurate. 
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